No enhancement of the localization length 
for two interacting particles in a random potential 
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We study two interacting particles in a random potential chain by means of the transfer matrix 
method. The dependence of the two-particle localization length A2 on disorder and interaction 
strength is investigated. Our results demonstrate that the recently proposed enhancement of A2 as 
compared to the results for single particles is entirely due to the finite size of the systems considered. 
This is shown for a Hubbard-like onsite interaction and also a long-range interaction. 
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In one dimension (ID), bosonization of the repulsive 
Hubbard model in the gapless phase shows that (onsite) 
disorder is a relevant perturbation and the ground state 
corresponds to a strongly localized phase B. Thus it 
came as a surprise, when Shepelyansky 0] recently ar- 
gued that the Hubbard interaction between two particles 
in a random potential would reduce the localization in 
comparison with independent particles. In particular, he 
obtained an enhancement of the two-particle localization 
length A2 independent of the statistics of the particles 
and of the sign of the interaction such that 
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in the band center. Here Ai is the single particle lo- 
calization length in ID and U the Hubbard interaction 
in units of the nearest-neighbor hopping strength. Shep- 
elyansky obtained this result by studying the matrix rep- 
resentation U of the Hubbard interaction in the disorder- 
diagonal basis of localized single particle eigenstates, i.e., 
U = UQ nmM = UY, jUi)<l>Ui)M*)Mi), with <f> n (i) 
the amplitude at site i of the single particle eigenstate 
with energy E n . Here, he assumed that <j> n (i) behaves as 
4>n{i) ~ Tn/yXi, where r n is a random number of order 
unity, and he neglected correlations among the r n at dif- 
ferent sites resulting in a Gaussian distribution for the 
matrix elements of Q. 

Support for this result was given shortly afterwards by 
Imry || with the help of a Thouless-type scaling argu- 
ment. Frahm et ai. ji] used the transfer matrix method 
(TMM) to study the two-particle Hubbard-Anderson 
problem without any approximations and have found nu- 
merically that A2 ~ A} 65 . They also measured the distri- 
bution of the matrix elements of Q and found a strongly 
non-Gaussian behavior. They attribute the deviation of 
the exponent from Eq. (Q) and the different distribution 
of matrix elements to the correlations among the r n ne- 
glected by Shepelyansky. Subsequently, an approximate 
calculation of A2 by Oppen et ai. j5| with Green func- 
tion methods lead to the hypothesis A 2 = Ai + C^-Af, 
with C ps 0.34 for bosons and 0.36 for fermions. Oppen 
et ai. also identified a scaling parameter U\\. Recently, 
Weinmann et ai. [|| argued that for a Gaussian matrix 



ensemble there exists a crossover from the \U\ behavior 
to the U 2 behavior for increasing U, but with a different 

exponent, such that A2 = Ai + Aj^^-X^ 2 for small U 
with W parametrizing the disorder. 

We feel that in this ongoing controversy about the in- 
teraction and disorder dependence, the system size has 
not been appropriately taken into consideration. We 
doubt the existence of the enhancement for an infinite 
system size, because there is only a vanishing probabil- 
ity that two localized particles are sufficiently close on 
an infinite chain to "see" each other. To substantiate 
this argument, we present in this Letter extensive com- 
putations of A 2 in dependence of the system size M. We 
first review the TMM approach to the two interacting 
particles (TIP) problem and suggest an improved TMM 
based on the self-averaging properties of the Lyapunov 
exponents. With this method we reproduce the A2 esti- 
mates of Ref. [Q for U = 1, system size M — 100, and 
all disorders. We find (cp. Fig. [I]), however, that the en- 
hancement A2/A1 decreases with increasing M. We also 
study the behavior of A2 for {7 = and show in Fig. |l| 
that it is not equal to Ai for finite M. But we recover 
the Ai value in the limit M — > 00. However, the en- 
hancement A2/A1 also vanishes completely in this limit 
as demonstrated in Fig. [j]. 

The Schrodinger equation for the TIP problem with 
Hubbard interaction is written in a suggestive form as 



lpn+l,m = [E - (e„ + e m ) - U5 n , m ] ip nt7n 



(2) 



where n,m = 1, . . . , M are the two site indices of the 
particles, E is the total energy of both particles, and ej 
is the random potential at site i. In the following, we use 
a box distribution [—W/2,W/2] for the ej. The single 
particle localization length in ID for such a distribution 
is known Q from second order perturbation theory in W 
to vary with disorder and energy as Ai » 24(4 — E 2 )/W 2 . 
The numerical results for E = based on the TMM 
lead to a slightly different prefactor 



Ai re 105/W 2 



(3) 



If one interprets (n, rn) as Cartesian coordinates on 
a finite lattice with M x M sites, the problem becomes 
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identical to a non-interacting Anderson model in 2D with 
disorder potential symmetric with respect to the diagonal 
n = m and with hard wall boundary conditions . One 
can rewrite Eq. (0) in the TMM form as 

with the symplectic transfer matrix 

r„=( E1 -*»- ff± "o 1 )' < 5 > 

where ip n = (V> n ,i> ■ ■ ■ , ip n ,m> ■■•> V'n.Af) is tne wave vec- 
tor of slice n, H± is the single-particle hopping term for 
the transverse (m) electron and {xn)i,m = (e« + fm + 
US njm )Si^ m codes the random potentials and the Hub- 
bard interaction [| . Note that in this approach the sym- 
metry of the wave function remains unspecified and we 
cannot distinguish between boson and fermion statistics. 

The evolution of the state is determined by the matrix 
product tjv = rin==i T n and we have 

(t + ')=™(&)- < 6 > 

Usually, the method is performed with a complete and 
orthonormal set of initial vectors (ipi,ipo) T . The eigen- 
values exp[27i(A^)] of (tLtn) 1 ^ exist for N — > oo due 
to Oseledec's theorem |1(J. The smallest Lyapunov ex- 
ponent 7 m in determines the slowest possible decay of the 
wave function and thus the largest localization length 
Amax = l/7min- We now define the two-particle local- 
ization length A2 as A max of the transfer matrix problem 

In the TMM studies of the Anderson Hamiltonian, one 
usually multiplies the transfer matrices in the n direction 
until convergence is achieved, e.g., for a 2D sample with 
strip width M = 100 and W — 2, this typically requires 
N = 5 x 10 6 multiplications for an accuracy of 1%. Thus 
one effectively studies a quasi- ID system of size M X N 
with M <C N. However, in the present problem, both 
directions are restricted to n,m < M. Iterating Eq. fl) 
only M times will not give convergence. Frahm et al. 0] 
have solved this problem in their TMM study by exploit- 
ing the hermiticity of the product matrix Qm = tJ^tm- 
Continuing the iteration (^) with t m , then with tm, and 
so on, until convergence is achieved, yields the eigen- 
values exp[2 7i (AT)] of (IT^Qm) 17 ^- This means that 
they repeat the same M X M sample N/M times. For 
M = 100 and W = 2 we found that N/M w 5000 is 
necessary to give an accuracy of 5%. Since this has to 
be done for every disorder configuration, the computa- 
tional effort is quite large. However, in the usual quasi- 
1D TMM for single particles, one may also halt the ma- 
trix multiplication every L -C N steps, compute an es- 
timate for 7(L) and finally average over all N/L such 
values [pr| . Due to the self-averaging property of the 



Lyapunov exponents |12j the averaged localization length 
equals the converged localization length X(N) within the 
accuracy of the calculation. This suggests the following 
method: We iterate Il^i 1 T r,M, where r enumerates dif- 
ferent samples of size M x M. Next we compute a local- 
ization length X(N) which is now already averaged over 
N/M different samples. This significantly reduces the 
computational effort and allows us to study much larger 
systems and achieve a significantly better accuracy than 
previously. 

In Fig. H we show the results obtained for the localiza- 
tion length of the TIP problem with Hubbard interaction 
by the present TMM. The data is calculated by configu- 
rationally averaging at least 10000 samples for each data 
point. This corresponds to an accuracy of 1% for the 
smallest W considered and better for larger W. We note 
that the data for U — 1 agree well with the results ob- 
tained in Ref. [| for 1.4 < W < 4 || and we find 
A 2 ~ W~ 3 8 . Furthermore, for W € [1.4, 5], the data can 
also be fitted reasonably well by 

A 2 = A x +A\f/(B + W), (7) 

with a = 2 with fit parameters A — 0.37 and B = 0.7. 
A fit with a = 3/2 according to Ref. || is considerably 
worse and without the denominator in Eq. (Q), as sug- 
gested in Ref. ||, the data cannot be fitted at all. We 
emphasize that the TIP data differ significantly from the 
2D TMM data with uncorrelated 2D random potential 

In Fig. H we have also plotted the behavior for {7 = 0, 
where the system reduces to two non-interacting parti- 
cles in ID. Before analyzing the data, let us first state 
our expectations: For {7 = 0, the two particles will 
localize independently at two arbitrary sites, say no, 
mo, with localization length Ai. The wave function 
can then be written as a product of two exponentially 
decaying single particle wave functions, i.e. ipn}m ~ 
exp[-|n - n \/Xi{Eo)} exp[-|m - m \/Xi(-E )}. Here, 
the eigenenergies are chosen such that E = 0. Two 
points are worth mentioning: (i) Let 71 (E a ) denote the 
inverse decay length measured by TMM for a single par- 
ticle with energy E a and 72(0) the inverse decay length 
for two particles with E = 0. For {7 = 0, the TMM 
for two particles will measure an inverse decay length 
72(0) ~ Y!ia 7 1 (Ea averaged over N pairs of states 
with energies E a and —E a . Since 71 (E 7^ 0) > 71 (0) as 
mentioned above Eq. (j|), we have 72(0) > 71 (0) and 
consequently A 2 (0) < Ai(0). (ii) The non-interacting 

two-particle wave function ipnj-n is not isotropic in the 
2D plane (n,m). Since the TMM will not necessarily 
measure the decay directly in the n direction, we expect 
that 72 will also contain information about the decay in 

other directions. These decay lengths of V>i°m are shorter 
and thus we again expect A2(0) < Ai(0). 

Figure || shows that contrary to our expectations the 
TMM does not under- but overestimate A2 as compared 
to Ai such that the enhancement X2(U = 1)/X 2 (U = 0) is 
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much smaller than the previously reported \ 2 (U = 1)/ Ai . 
Also included in Fig. g is a ID TMM result and we see 
that for W > 6 — where the perturbative result (||) is no 
longer accurate — both the U = and the U — 1 data 
agree quite well with Ai. For W G [1.4, 4] a power law fit 
gives \ 2 (U = 0) ~ W 3 - 6 . A fit with a = 2 as in Eq. (@) 
can also describe the data for small W with A = 0.19 
and B = 0.5. 

We emphasize that the deviation from Eq. (g) of our 
data for vanishing U does not result from an effective cou- 
pling which was suspected [Q to be introduced into the 
problem because of numerical instabilities of the TMM. 
Our tests did not show such a coupling. We also note 
that the two-particle TMM data is quite different from 
data of a 2D single-particle calculation for M = 100 as 
shown in Fig. |^. 

Studying the U dependence of A2, we note that 
X 2 (U) = Xz{—U) at E = due to particle-hole sym- 
metry. Consequently, we only show data for U > in 
Fig. H We observe neither the proposed U 2 nor the \U\ 
behavior. Rather the data obey \ 2 {U) - A 2 (0) w U 63 
for small U. Also, A2(0) ~ 8.45 does not agree with 
Ai ~ 105/16, as discussed above, and we cannot identify 
a scaling parameter as suggested in Ref. . 

From Fig. || we see that for W > 6, the data for [7 = 
and U = 1 coincide with the single-particle TMM data 
already quite well. The enhancement A2/A1 is large only 
for small disorder (W < 4). However, in this region, 
the computed values of A2 become comparable to the 
size of the system such that an increasing part of the 
wave function is reflected at the hard wall boundaries. 
This leads to an overestimation of X 2 . In order to dis- 
tinguish this artificial from the proposed enhancement 
due to interaction alone, we have therefore systematically 
studied the finite-size behavior of A2 at various disor- 
der values and interaction strengths. As an example, we 
consider the case W — 3. Here the localization lengths 
for X 2 (U = 0) « 19.1 and \ 2 (U = 1) w 25.5 may just 
be small enough for the wave function to fit into the 
100 x 100 box without too strong boundary effects. Fur- 
thermore, the enhancement is still sufficiently large to be 
numerically detectable. In Fig. [I], we have shown A2 as a 
function of M ~ x l 2 . Our reason for plotting the data as 
a function of M -1 / 2 is pragmatic: a plot versus M _1 or 
TV/- 1 / 4 shows a downward or upward curvature, respec- 
tively, for the largest M. The curvature is smallest for 
A/ -1 / 2 and the extrapolation towards infinite M is most 
reliable. We see in Fig. |l| that (i) the data for U = 
approach the single-particle result Ai « 105/9 for large 
AI, and, more significantly, (ii) the data for U = 1 also 
approach this non-interacting ID result. Therefore the 
overestimation of X 2 (U = 0) vanishes for large M, sup- 
porting the validity of the TMM approach. However, in 
this limit the enhancement X 2 (U = 1)/Ai vanishes com- 
pletely, too. We note that the observed decrease in the 
localization length A2 is quite different from a 2D TMM, 
where the localization length is known to increase with 
increasing system size |J. We attribute this difference 



to the symmetry of the TIP disorder potential which is 
equivalent to an effective long-range correlation which 
furthermore increases with increasing M. 

We have also studied the behavior of the two-particle 
localization length for a system with long-range interac- 
tion, since Shepelyansky argues in Ref. (ll| that long- 
range interaction should also lead to the enhancement of 
A2. The Schrodinger equation is given as 

4>n+i, m = [E- (e„ + e m ) - U/(\n -m\ + 1)] V>n,m 

— VVi,m+l — 1pn,m-l ~ fn-Lm, (8) 

such that the on-site interaction is equal to the Hubbard 
interaction. Clearly such a system should exhibit even 
stronger boundary effects and we thus expect an even 
stronger enhancement of A2 for finite systems. In Fig. ^, 
we have also shown A2 for the long-range interaction with 
U = 1 at E = and the same accuracy of 1% as for the 
Hubbard interacting system. As expected, we observe 
that the enhancement for all M is stronger than for the 
Hubbard interacting model. But this enhancement again 
vanishes for M — > 00. 

In summary, we have studied the interaction-induced 
enhancement of the localization length for two particles 
in a ID random potential by means of an improved TMM. 
For U = 1 and M — 100, the results agree well with 
previously published data. However, the U — data 
is quite surprising since only in the limit M — > 00 do 
we recover the expected non-interacting ID result Ai. 
However, in this limit the enhancement for finite U also 
vanishes. Therefore we must conclude that the transfer 
matrix method applied to the ID Anderson model for 
two interacting particles measures an enhancement of the 
localization length which is entirely due to the finiteness 
of the systems considered. This enhancement might be 
relevant for applications in mesoscopic systems. 
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FIG. 1. Two-particle localization length A2 as a func- 
tion of the system size at W — 3 for Hubbard (□) and 
long-range interacting systems (*). Dashed lines are extrap- 
olations from the data for M £ [80, 170]. The 5 leftmost 
data points for U — and Hubbard U = 1 corresponding 
to M = 180, 200, 240, 300, 360 have been computed with only 
5% accuracy and have not been considereded for the extrap- 
olations. The diamond (o) indicates the value of Ai. 
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FIG. 2. Two-particle localization length A2 (accuracy at 
least 1%) at energy E = for system size M = 100 and 
interaction strength U = (o) and 1 (□). The solid line is 
given by Eq. (^). Single-particle TMM data with 1% accuracy 
for ID chains (x) and 2D strips of width M = 100 (+) are 
also shown. 
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FIG. 3. U dependence of A2 for W — 4. The dashed line 
represents the power law fit \2(U) — A2(0) ~ U 63 . The lower 
curve shows the same data on an enlarged U scale. Arrows 
indicate the corresponding axes. 
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